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We show that it is possible for chaotic systems to display the main features of coherence resonance. In 
particular, we show that a Chua model, operating in a chaotic regime and in the presence of noise, can exhibit 
oscillations whose regularity is optimal for some intermediate value of the noise intensity. We find that the 
power spectrum of the signal develops a peak at finite frequency at intermediate values of the noise. These are 
all signatures of coherence resonance. We also experimentally study a Chua circuit and corroborate the above 
simulation results. Finally, we analyze a simple model composed of two separate limit cycles which still exhibits 
coherence resonance, and show that its behavior is qualitatively similar to that of the chaotic Chua system. 



When a dynamical system is subjected to an external pe- 
riodic forcing, it is a standard result that synchronization be- 
tween the system and the forcing can occur under a large va- 
riety of conditions. A resonance is defined as the presence 
of a maximum in the response of the system as a function 
of some control parameter (for instance, the frequency of the 
external signal). Although one might naively believe that fluc- 
tuations, either in the forcing or in the intrinsic dynamics, will 
worsen the quality of the synchronization, it is nowadays well 
established that, in some cases, the response of a nonlinear 
dynamical system to an external forcing can be enhanced by 
the presence of noise (fluctuations). The prototypical and pio- 
neering example is that of stochastic resonance by which 
a bistable system under the influence of a periodic forcing, and 
in the presence of fluctuations, shows an optimum response, 
a resonance, for a given value of the noise intensity. The rel- 
evance of this phenomenon has been shown for some physi- 
cal and biological systems described by a nonlinear dynamics 



That noise can have a constructive role has been one of 
the most astonishing discoveries of the last decades in the 
field of stochastic processes. Besides the above mentioned 
stochastic resonance, purely temporal dynamical systems can 
display phenomena such as noise-induced transitions [||] or 
noise-induced transport [^. In spatially extended systems, on 
the other hand, noise is known to induce a large variety or or- 
dering effects [|^]. In all the cases, the common feature is that 
some sort of order appears only in the presence of the right 
amount of noise. 

The possibility of having stochastic resonance without the 
need of an external forcing has attracted much attention re- 
cently [p[-pl]|. In particular, the phenomenon named coher- 
ence resonance Jill ] was shown to appear in excitable sys- 
tems under the influence of fluctuations. An excitable sys- 
tem has a stable fixed point with a finite basin of attraction. 
When a perturbation is such that the system crosses a thresh- 
old value, the return to the fixed point is by executing a large 
excursion in the configuration space, thus generating a pulse 



in the time evolution. One of the main features of excitable 
systems is that the generated pulse is basically independent 
of the magnitude of the perturbation that induced its firing. 
Therefore, the duration of the pulse, the excursion time t^, 
is a characteristic of the system and not of the perturbation. 
The total time between pulses, tp, is composed of two times: 
the excursion time te and the time needed for the activation 
of the pulse, ta- If the firing of the pulses is produced by 
random fluctuations, the activation time ta is a random vari- 
able. According to Kramers formula, for small noise, the 
mean activation time behaves as (to) ^ cxp{A/ D^) where 
yl is a constant and D is the noise intensity [[T^. The vari- 
ance of the activation time is (T^[ta] ~ (^a)^- At the same 
time, the excursion time depends weakly on D, such that 
its mean value (te) can be considered constant and its vari- 
ance can be estimated as cr^[te] ~ D'^{te)- For small D, we 
have (ta) ^ (ie) and we can approximate the time between 
pulses by the activation time tp w ta- The relative fluctua- 
tions of the time between pulses, defined as R — a[tp]/{tp), 
is in this limit of small noise R w CT[ia]/ {ta) ~ 1- For large 
noise, the activation time is very small and the system fires 
a pulse every time it returns from an excursion. The pulse 
time is dominated by the excursion time and we can approx- 
imate R « a[te]/{te) ^ D{te)^^^^. If the excursion time 
is large, and the threshold of excitation is small, it is possible 
that for intermediate values of D, it is R{D) < 1. In this 
case, and according to the generic behaviors described above 
(R{D) -> 1 for small D and i?(L>) - ^(ie)-!/^ for large £>) 
there will be a minimum in the relative fluctuations of the time 
between pulses. This is the signature of coherence resonance. 
A similar effect is that of stochastic resonance without exter- 
nal periodic force which can occur in a system near a limit 
cycle bifurcation point 10 1. 



To summarize, the main feature of a system displaying co- 
herence resonance is that a quasi-periodic signal is generated 
by a combination of the internal nonlinear dynamics and fluc- 
tuations without the need for the presence of an external, de- 
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terministic, periodic signal. The periodicity of the pattern is 
optimal (resonance) for a certain value of the noise intensity. 
The original studies have been extended to consider other ex- 
citable systems such as the FitzHugh-Nagumo model [|l3|jl4|, 
the Hodgkin-Huxley model for neurons [ p3[ ] and the Yamada 
model for a self-pulsating semiconductor laser jl^]. Coher- 
ence resonance has been also observed in dynamical systems 
close to the onset of a bifurcation [17| as well as in other 
bistable and oscillatory systems jr^,19]. Experimental evi- 
dence for the existence of coherence resonance has been given 
for a laser system and for excitable electronic circuits 

In this paper we will prove that it is possible to display 
the main features of coherence resonance in chaotic and other 
bistable systems in which the attractors are not of the fixed 
point type. Although we believe that our results are quite gen- 
eral, we will consider specifically a Chua circuit operating in 
a chaotic regime with two independent, symmetric, attractors. 
We will argue that the existence of a well defined character- 
istic time when moving around each attractor is a necessary 
ingredient for the ocurrence of coherence resonance. This 
characteristic time plays the role of the excursion time for ex- 
citable systems in the sense that only during a small fraction 
of this time, when the trajectory comes as close as possible to 
the other attractor, the fluctuations can induce jumps between 
the two attractors. We will give numerical and experimental 
evidence that such a Chua circuit in the presence of noise, can 
undergo oscillations whose regularity is optimal (in a sense 
to be precisely defined later) for some intermediate value of 
the noise intensity. Later, we will also show that the basic 
ingredients for this new kind of coherence resonance are al- 
ready present in a simpler toy model with two separated limit 
cycles. That the addition of noise can induce some degree 
of regularity in a chaotic system has been shown recently in 
a different context related to the existence of noise-induced 
sychronization of chaotic systems [ p3[ ] 

Let us consider the Chua system, in its dimensionless form, 
under the presence of additive noise [e3]; 
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FIG. 1 . Time series for the x variable of the Chua system given by 
Eqs. (§) for three different noise levels: a) D = 0.02, h) D = 0.08, 
optimum noise level, and c) D = 0.16. 
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where ^(t) is Gaussian white noise, of zero mean and correla- 
tions {i{t)^{t')) = D'^S{t - t'). The nonUnear function h{x) 
is given by h{x) = bx + ^^{\x + 1| — |x — 1|). We have 
taken the values a = -1/7, = 2/7, a = 4.60, /? = 6.02, 
7 = 0, for which the Chua system has two chaotic attractors: 
a single scroll and its mirror image. Depending on the initial 
conditions, the system will rotate around one attractor or the 
other. In other words, in the absence of fluctuations, the attrac- 
tors are independent and trajectories can not jump from one to 
the other The movement around each attractor has a well de- 
fined mean angular frequency wq, which for these values of 
the parameters is « 3. 
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FIG. 2. Standard deviation normalized by the mean time u / (T) 
for the Chua system given by Eqs. (^ for noise levels ranging from 
0.01 to 0.17. 

In figure ^we plot three trajectories of the variable x(t) cor- 
responding to increasing levels of the noise intensity D. We 
observe three qualitatively different behaviors when increas- 
ing the noise level. When D is very small (figure |l]a) the aver- 
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age residence time, i.e. the time between jumps, is large and 
the system spends most of the time rotating around one of the 
attractors. For this situation, the dispersion of the residence 
time is also large. As D increases, approaching an optimum 
value (figure |l]b), the system jumps between the two attrac- 
tors more regularly. These jumps occur, as already mentioned, 
approximately when the trajectory passes closest to the other 
attractor. Finally, when D is very large, the system jumps 
more often but these jumps may start from different points on 
the trajectory and the behavior of the system is more irregular 
(figure The enhanced regularity that occurs for interme- 
diate values of the noise can be clearly observed in figure ^ 
where we plot the standard deviation, cr\T], of the residence 
time in each attractor, normalized to its mean value, (T), as a 
function of the noise intensity. This curve exhibits a minimum 
at a noise level D « 0.08. The presence of this minimum is 
the clearest signature of coherence resonance. Another ob- 
served indicator of coherence resonance is the existence of a 
peak in the the power spectrum S{f ) of the signal at a finite 
frequency [ n , 1^ . Moreover, our data show a maximum in 
the ratio between the height and the width of the peak of S{f) 
at Z? ~ 0.07. Similarly, we observe that the time-correlation 
function, C(t), has the longest tail and the lowest minimum at 
values of the noise close to the optimal level. 

Although there is a good correspondence between a real 
Chua circuit and the system of differential equations (|l]), it 
is obvious that the numerical results deal necessarily with an 
idealized Chua circuit. In order to analyze the robustness of 
the observed phenomenon, we have performed experiments 
in a real Chua circuit constructed according the classical de- 
sign (see e.g. ref [p4[]) with the following parameters |p5[]: 
Ci = 20 nF, C2 = 100 nF, R = 1100 Vl, a ^ -1/7, b = 2/7. 
The noise has been generated with a standard Hewlett-Packard 
function generator and its intensity has been varied from zero 
to a few volts. As in the numerical study, the original time 
series for the x{t) variable has been converted into a variable 
u{t) taking the values +1 and —1 for each of the attractors. 
Using this variable we have computed the normalized stan- 
dard deviation of the residence time, cr [T] / (T) . In figure |^ we 
plot this quantity as a function of the noise intensity. Again, a 
clear minimum appears for an optimal value D k, 5500 mV, 
so confirming the numerical results. In figure ^ we plot the 
power spectrum of the digital variable u{t). We notice the de- 
velopment of a peak at a finite frequency for an intermediate 
noise level giving further evidence of a regular behavior. 

To gain more insight into the dynamics of this chaotic sys- 
tem we consider a simplified two variable system {xi^X2) 
with two stable limit cycles. The first one, Ci, is around the 
unstable fixed point (1,0) and the second one, C2, around the 
unstable fixed point (—1, 0). There are no other stable fixed 
points or limit cycles in the system. We assume that the limit 
cycles are circumferences of radius R close to but smaller than 
1, and that they have a constant angular speed tJo- Under these 
circumstances, which limit cycle is chosen as a dynamical at- 
tractor depends exclusively on the initial condition p^ ] . Let 
us add now some noise to the dynamics. If the noise inten- 



sity is small, the modification to the trajectories will be small. 
Moreover, the probability that noise induces a jump between 
the attractors is only significant near the closest points in the 
limit cycles, i.e., when the trajectory passes closest to the ori- 
gin of the coordinate system. If the system does not jump at 
this point then it has to wait for a complete rotation for an- 
other chance to jump. Hence, the rotation period 2-k jtOQ plays 
the role of the excursion time in the excitable system in 
the sense that, for moderate levels of noise, the system cannot 
jump to the other attractor during this time. 
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FIG. 3. Standard deviation normalized by the mean time a / (T) in 
the case of the Chua circuit (see the text for details of the parameters) 
for noise levels ranging from 2000 to 7500 mV. 
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FIG. 4. Power spectrum S{f) of the digital signal u{t) for the 
Chua circuit for different noise levels: D — 2000 mV (squares), 
D = 4000 mV (rhombi), D = 6000 mV (triangles). The presence 
of a peak at a finite frequency is an evidence of coherence resonance. 

Let us define the variable u — sign(a;i). This is equal to +1 
when the system is in attractor Ci and — 1 when in attractor 
C2. The evolution of u{t) can be described by a series of time 
intervals alternating the values of +1 and —1. According to 
the previous argument, for small noise, the duration T of each 
time interval is a random variable taking values which are an 
integer multiple of 2tt/ujo. The probability that the system 
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jumps between the attractors exactly after n cycles follows 
the geometric distribution: P{T = 2TTn/u!o) = p{l — 
where p is the probability that the system jumps between the 
two attractors during the time they are closest at each cycle. 
This probability p will be a function of the noise intensity D 
and also of the angular speed ujq determining the time the sys- 
tem is ready for jumping. Using this distribution, the relative 
fluctuation of the time betweenjumps is (j[T]/ (T) = ^/l — p. 
Therefore, for small D, p will be small and R will initially 
decrease with D. Since, according to the general argument of 
[ pi] ] developed at the beginning, R will eventually grow for 
large D, we expect a minimum in a plot of (t[T]/{T) versus 
D. When reexamining figures^ and g, it can also be observed 
that both limits, low and high noise intensity, behave roughly 
as estimated for the simpler model described above indicat- 
ing that our simplified model captures the main ingredients of 
coherence resonance in the more complicated chaotic system. 
However, there is an important difference between the simple 
dynamical model and the Chua system. In the former, the tra- 
jectory is always on the limit cycle and the system is ready to 
jump to the other attractor at any cycle. On the contrary, in 
the Chua system, after the trajectory jumps from one attractor 
to the other, the motion usually starts close to the center of the 
attractor, in an inner orbit, and it is not ready to jump to the 
other attractor until the outer orbits are reached. This fact can 
be guessed from figure |l] when looking at the times the system 
jumps from one attractor to the other This difference might 
be responsible from the fact that in the Chua model R does not 
clearly approach to 1 as Z? is decreased. On the other hand, 
if noise is arbitrarily increased, the Chua system saturates to a 
unique limit cycle thus losing its characteristic behavior. 

In conclusion, we have shown, both numerically and ex- 
perimentally, that coherence resonance can be observed in a 
chaotic system. We have also shown that a simple model, 
composed of two separate limit cycles, is able to exhibit co- 
herence resonance. Within this model we were able to predict, 
for instance, the limits of the normalized standard deviation of 
the residence time by a simple analytical approximation. The 
behavior of the chaotic Chua system follows qualitatively the 
results derived in the simple model, with coherence resonance 
illustrated by the dependence of several different quantities on 
the noise intensity. Finally, we consider particularly interest- 
ing the fact that the combination of noise and chaos can lead 
to some degree of regularity in the system. 
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